Detonative Propagation and Accelerative Expansion of the Crab Nebula Shock Front 
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The accelerative expansion of the Crab nebula's outer envelope is a mystery in dynamics as 
a conventional expanding blast wave decelerates when bumping into the surrounding interstellar 
medium. Here we show that the strong relativistic pulsar wind bumping into its surrounding nebula 
induces energy-generating processes and initiates a detonation wave that propagates outward to form 
the current outer edge, namely the shock front, of the nebula. The resulting detonation wave, with 
a reactive downstream, then provides the needed power to maintain propagation of the shock front. 
Furthermore, relaxation of the curvature-induced reduction of the propagation velocity from the 
initial state of formation to the asymptotic, planar state of Chapman- Jouguet propagation explains 
the observed accelerative expansion. The essential role of detonative propagation in the structure 
and dynamics of the Crab nebula offers potential richness in incorporating reactive fronts in the 
description of various astronomical phenomena. 

PACS numbers: 95.30.Lz, 98.58.Mj, 47.40.Rs 



The Crab nebula is the most observed supernova rem- 
nant since its birth in A.D. 1054. Its explicit birth date, 
together with accurate observations of its subsequent 
evolvement through optical [ll-Q and radio 0] facilities 
in the last several decades, has yielded rich information 
on its structure and dynamics. In particular, it is shown 
that this nebula has two major features, namely opti- 
cal line-emitting filaments, which are supposed to be the 
dense condensations of the thermal gas @, and an 
expanding radio synchrotron bubble shell, which is rec- 
ognized as the outer edge of the nebula composed of light 
relativistic gas and magnetic fields Furthermore, us- 
ing observations from 1939 to 1992, and by assuming a 
constant expansion speed for the nebula, the date of the 
supernova outburst was estimated to be A.D. 1130±16 
Q, which is 76 years after the historically observed date 
of A.D. 1054. This discrepancy then implies that the 
nebula expansion is actually accelerative. A similar re- 
sult was obtained through an analysis of the nebula's 
radio shell Q, at « 2 pc = 6 X 10 18 cm. Based on the 
most recent result that the current expansion velocity of 
the nebula shell is v p ~ 1.5 x 10 8 cm s _1 , and if a con- 
stant acceleration v c is assumed, one can then estimate 
that the acceleration is 0.82 x 10 -3 cm s~ 2 . 

The above observation, however, cannot be satisfacto- 
rily explained by current understanding of the dynamics 
of the expanding front, on two counts, both based on the 
energy source needed to sustain the front propagation. 
First, it has been suggested that the energy from the spin 
down of the central pulsar provides the needed energy 
source 0, @]. However, it is not clear as how the spin- 
down energy of the central pulsar, carried by the pulsar 
wind, can be transferred to the kinetic energy needed to 
support the expansion of the nebula's outer edge. This is 
because the matching radius of the strong MHD shock, 



r s w 3 x 10 cm = 0.1 pc, formed by the highly rela- 
tivistic pulsar wind bumping into the surrounding neb- 
ula, is much smaller than the outer edge of the nebula, 
r<j w 2 pc, to induce effective coupling (see Fig. 1 for an 
illustration of these radii). Second, present theoretical 
considerations generally describe the outer edge of the 
Crab nebula as the front of a (magneto) hydrodynamic 
blast wave that bumps into the surrounding interstellar 
medium (ISM) [7]. However, since a conventional ex- 
panding blast wave decelerates Q as it transfers kinetic 
energy to its surrounding, such a feature not only contra- 
dicts the accelerative nature of the nebula's outer shell, it 
cannot even support the steady expansion of the front. It 
is emphasized that the acceleration, while low in magni- 
tude, being just ~ 5% in the entire life of ^1,000 years, is 
nevertheless finite and as such accentuates the observed 
absence of deceleration. 

In view of the above concerns, it is of intrinsic inter- 
est to explore whether there are other energy sources and 
coupling mechanisms for the observed accelerative expan- 
sion of the nebula shock front. We suggest herein that 
the shock front is actually a detonation wave, which is 
a shock wave with an exothermic, reactive downstream. 
The exothermic heat release then provides the needed 
energy for the sustenance of the propagation. Further- 
more, relaxation of the front curvature upon expansion 
facilitates the front propagation, thereby resulting in the 
observed accelerative expansion. We now present our 
interpretation and implementation of the detonative de- 
scription of the nebula shock front propagation, with cor- 
responding estimates of the various dynamic and thermo- 
dynamic parameters supporting such a mechanism. 

Detonation wave at the nebula envelope To facilitate 
appreciation of the role of the detonation wave on the 
problem at hand, it is instructive to first discuss its struc- 
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FIG. 1: Schematic structure of the Crab nebula. The cross 
in the center represents the pulsar. The region immediately 
around the pulsar is the pulsar wind domain, bounded by the 
surrounding nebula with an interface of strong MHD shock 
at r s =0.1 pc. Beyond this radius is the nebula region, with 
an outer edge at rd = 2 pc, where the detonation wave forms 
as the nebula expands into the ISM. 

ture, initiation and propagation based on understanding 
largely gained from studies of reactive gasdynamics d,Q. 
Briefly, in the wave stationary frame, a detonation wave 
consists of a leading hydrodynamic shock which slows 
the incoming supersonic flow to subsonic and changes its 
thermodynamic states according to the equation of state 
of the medium. This leads to the Neumann state imme- 
diately downstream of the shock. For conventional mate- 
rials including that of interest in the nebula, the density, 
pressure and temperature are increased upon crossing the 
shock, which induces reaction of the material after an ig- 
nition delay length, This forms the ZND structure 
of the detonation wave. For reactions that are highly 
temperature sensitive, the subsequent emergence of the 
reactions is abrupt and basically complete, in the man- 
ner of a step function. Figure 2 shows a schematic of the 
profiles of the various properties. 

Propagation of a detonation wave is affected by the 
downstream rarefaction waves that continuously over- 
take the shock front and weaken its strength until the 
velocity downstream of l ig becomes sonic. For ID pla- 
nar propagation without loss, the sonic state is attained 
at downstream infinity and steady-state propagation is 
attained, resulting in the Chapman- Jouguet (C-J) mode 
of propagation. This feature of detonation propagation 
therefor assures the non-deceleration of the Crab nebula 
shock front. The acceleration, however, is due to the 
relaxation of the so-called curvature effect. Specifically 
for an expanding front the divergence of the streamline 
downstream of the shock slows down the flow velocity, 
advances the sonic state from far downstream towards 
the front, and consequently reduces the propagation ve- 
locity from the C-J value, with the reduction increasing 



with increasing front curvature. This curvature-induced 
reduction has three implications. First, there exists a 
minimum, critical radius R cr for the initiation of the deto- 
nation propagation. Second, once initiated and the front 
curvature relaxes as it expands, the propagation velocity 
would continuously increase to approach the C-J value. 
Third, a minimum, critical energy is needed for a deto- 
nation to attain its critical radius. 

We shall now apply the above understanding to ana- 
lyze the propagation of the nebula front. Referring to 
Fig. 1, it is proposed that from the beginning of the 
nebula's evolution, the pulsar wind has been pushing its 
surrounding medium outward continuously until the in- 
terface expanded to r s , at which a forward and a re- 
verse detonation fronts are simultaneously formed, with 
the forward front propagating outward and the reverse 
front propagating inward. The continuously expanding, 
forward detonation wave then forms the outer envelope 
of the Crab nebula observed today. The motion of the 
latter, however, is rapidly arrested by the outwardly di- 
rected pulsar wind, forming the MHD shock. The radius 
r s , at which the reverse detonation front is stabilized, can 
thus be identified as the critical radius for detonation ini- 
tiation, which we shall discuss in following context. 

Our analysis is based on a spherically symmetric model 
describing the roughly spherical nebula's outer shell as 
was observed As a first attempt to the detonation 
interpretation of the nebula's expanding shell, we do not 
consider the filament structure of the nebula region 0, Q , 
but will make some observations of its possible effects 
later. We further note that previous analysis Q shows 
that in order to convert the kinetic energy flux effectively 
into the observed synchrotron luminosity, the ratio of 
the magnetic field energy to the particle kinetic energy 
should be very small. Then it is reasonable to neglect the 
magnetic field when considering the energy components 
in this wind- nebula shock at r s , as well as in the more 
distant nebula regions where the detonation takes place. 
In addition, we also assume that the gravitational energy 
is much smaller than the kinetic energy of the flow. 

With the above assumptions, the problem simplifies 
to that of the structure and propagation of a spheri- 
cally expandin g de tonation wave which have already been 
well analyzed (ToL Ulj ]. and was presented in [&]. We 
shall therefore simply state the key results relevant for 
the present problem. Through estimates of the various 
dynamic and thermodynamic parameters governing the 
present process, we shall demonstrate that the observed 
phenomena associated with the expanding nebula front 
are consistent with the description based on detonation 
theory. 

Parametric structure of the detonation wave We de- 
fine the radius R s , the propagation speed D, the heat 
of reaction q c and the reaction induction length ^ g for 
the detonation wave. The corresponding dimensionless 
variables are k g = ?i g /i? s , q c = g c 7/a^ and D = D/a u 
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the Mach number of the detonation velocity, where a = 
\pjpfp is the sound speed, 7 is the adiabatic index and 
the subscript u denotes the upstream ISM. Focusing on 
the evolution of the detonation Mach number, we can 
express it as D — -Dcj(1 — A)> where A <C 1 is a small 
number representing its deviation from the asymptotic 
C-J value caused by the curvature stretch. Then we have 



.9] 



Dh = 2( 7 + 1) + 2( 7 2 - 1)^ « 2( 7 2 - 1)^, (1) 



8 7 2 



r 



(2) 



for large values of q c relevant for the present situation of 
strong detonation, which has an estimated Mach number 
of 170 (see Table I). The typical temperature of the un- 
burned upstream ISM is [l2j T u ~ 10 4 K and the density 
is p u ~ 0.5 cm~ 3 . Other parameters of the upstream ISM 
are readily calculated and listed in Table I. Furthermore, 
according to eq. (1), the heat of reaction that supports 
the expansion of the detonation wave is estimated as 



5.9 x 10 i2 J/mol = 61 MeV/atom, 



(3) 



which is relatively large compared to typical values of 
light nuclear reactions and those of chemical reactions. 

It is also noted that in the above calculations, the adi- 
abatic index of the upstream ISM gas has been assumed 
to be 7 = 1.1, resulting in a large thermal capacity, 
c p « (v « 10i?, with R being the molar gas constant. 
The index 7 = 1.1 is intermediate of isothermal (7 = 1) 
and adiabatic (7 = 5/3 for monatomic gas) fluids, and is 
a typical value for ISM with partially effective thermal 
conduction. Furthermore, detonation cannot occur for 
isothermal ISM gas with 7 = 1, while for the adiabatic 
limit case of 7 = 5/3, the reaction energy release q c will 
decrease to ~l/6 its current value, and the velocity de- 
viation A will decrease to 1/4 its current value. These 
variations will cause tuning of the detonation parame- 
ters, but will not change the principal conclusions made 
in this letter. 

We can further calculate the physical parameters de- 
scribing the Neumann state just downstream of the shock 
front but before the reactions are initiated After the 
reactions have taken place, the downstream (burned) gas 
will eventually approach the C-J state, whose parame- 
ters can also be readily determined. Parameters of the 
upstream state, the Neumann state and the burned C-J 
state are listed in Table I and schematically shown in Fig. 
2. 

It is noted that under the Neumann temperature Tn = 
1.4 x 10 7 K, some nuclear reactions, e.g. the p-p cycle 
fl3j | . are able to be initiated. However, as the density in 
the nebula region is smaller by a factor of ~ 10 -25 com- 
pared to that in the stellar central regions, reaction rates 



TABLE I: Gas parameters of the detonation front: upstream 
ISM (u), Neumann state (N) and burned C-J state (C-J). 
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FIG. 2: Temperature, density, velocity and pressure of the up- 
stream ISM, the compressed Neumann state and the burned 
C-J gas are shown as step functions in each panels. The Neu- 
mann state crosses the induction length of h g = 1.0 x 10 -4 pc, 
while the upstream and burned gas extends outward to infi- 
nite and inward to the inner edge of the nebula, respectively. 



of these nuclear processes are very small. Since typical 
heats of chemical reaction are also very small, alterna- 
tive energy-generating processes need to be identified to 
supply the energy of q c * 61 MeV/atom. One possible 
energy source for the detonation is the inverse Compton 
process [Tij . which is a commonly used mechanism 
in astrophysics. Crab nebula's high energy 7-ray emis- 
sions and flares are supposed to be produced by inverse 
Compton scatterings and very probably, part of 

the high energy photons will interact with the nebula gas 
to support expansion of the envelope. 

Initiation of the detonation wave The standing in- 
ward shock front (see Fig. 1 and reference Q) is assumed 
to be the radius where the forward and reverse detona- 
tions began to form and separated in the early days of 
the nebula evolution. Since the reverse shock is expected 
to be immediately captured by the pulsar wind, the ra- 
dius of its front varies little since the formation of the 
binary detonations. Consequently the critical radius for 
the formation of the forward propagating detonation can 
be identified as that of the standing MHD reverse shock, 
namely R CI = r s =0.1 pc. 

Directly referring to eq. (14.8.29) in |9| that relates the 
induction length and the critical radius, and assuming 
(0 = 4, which is consistent with the reaction activation 
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temperature estimated based on the Neumann state, the 
induction length can be estimated as 



R c 



1.0 x 10~ 4 pc, 



(4) 



which is a small number compared to the radius of the 
detonation wave, = 2 pc, and is therefore consistent 
with the requirement of the formulation. We then have 
the quantitative expression of eq. (2): 



8 7 2 



r 



i 



0.004 pc 



(5) 



Attainment of the critical radius for the detonation 
initiation requires the initial blast wave to possess a 
corresponding critical energy E c 
(14.8.30) in [g, we find: 



181 ] . Referring to eq. 



E r 



7 2 _ i 



3+1 



3 x 10 45 erg, (6) 



where j ' — 2 and fcj = 5.31 for the spherical geometry. 
This critical energy is much smaller than the part of pul- 
sar spin-down energy carried by the pulsar wind, esti- 
mated to be ~ 10 49 erg 0, UK]. So the energy carried 
by the relativistic pulsar wind is sufficient to initiate the 
forward detonation wave. For those supernova remnant 
systems whose pulsar wind energies are not large enough 
to reach the critical value of E CI ~ 3 x 10 45 erg, detona- 
tion cannot be successfully initiated; instead, it will de- 
celerate as it defuses into the surrounding medium. We 
can therefore classify those nebulae around supernovae 
into two groups: namely acceleratively expanding ones 
led by detonation waves, and defused ones formed by de- 
celerating blast waves. 

Accelerative expansion of the detonation front Subse- 
quent to attaining the critical radius, the forward deto- 
nation then propagates with the velocity 

D = £> C j(l-A) = 1.5xl0 3 x(l- °-° 04 P ° )km s' 1 (7) 

in the nebula domain that begins at r s = 0.1 pc and 
ends at the current nebula envelope at rd = 2 pc. The 
variation of the nebula envelope velocity as a function of 
the envelope radius is shown in Fig. 3, demonstrating 
an accelerative expansion of the wave front. It is seen 
that at the early age of the nebula, when R s is slightly 
larger than 0.1 pc, the propagation velocity of the deto- 
nation is smaller than its present value by about 5 per- 
cent. Furthermore, the detonation velocity accelerates 
more rapidly in the early days of the front formation, 
and gradually approaches the C-J velocity as it evolves. 
We have therefore successfully explained the accelerative 
propagation of the nebula front. 




FIG. 3: Increase of the detonation propagation velocity as a 
function of the wave front radius. The detonation velocity is 
~ 1425 km s _1 at R B = 0.1 pc where the detonation begins 
and gradually approaches the C-J velocity of ~ 1500 km s _1 
as the wave expands. 



Summary and Discussion In this first attempt to in- 
tegrate detonation features in the structure and dynamics 
of the Crab nebula, we have demonstrated that the con- 
cept is a viable one in terms of the observed phenomena 
as well as the various estimated dynamic and thermo- 
dynamic parametric values associated with the descrip- 
tion. In particular, the following understanding has been 
gained. (1) There is a supersonic detonation wave at 
the front of the expanding nebula with a Mach number 
of ^170; the wave is supported by downstream energy- 
generating reactions. (2) The surrounding ISM is com- 
pressed to a temperature of 10 7 — 10 s K and a density 
of 1 — 10 cm~ 3 by the initial shock and the subsequent 
reactive heating. (3) The envelope of the nebula is ac- 
celerating due to relaxation of the front curvature as it 
evolves to the planar limit upon continuous expansion. 

We close this discourse with some additional consid- 
erations. First, there is another possible mechanism in- 
volving the pulsar wind energy that may account for the 
accelerative expansion of the Crab nebula outer envelope. 
That is, the downstream reaction of the inward detona- 
tion may last for a long time, in which case the delayed 
exothermic reaction takes place close to the outer enve- 
lope and provides the kinetic energy needed for its ac- 
celerative expansion. This mechanism is similar to the 
flame acceleration in obstructed channels [20} . However, 
it is an open issue to confirm the existence of delayed 
exothermic reactions in the Crab nebula. 

Second, it is noted that the filament structure of the 
nebula could be related to the well-established cellular 
morphology of detonation fronts Q . These cells are con- 
sequences of interacting triple-shock ensembles, which 
in turn generate complex local subsonic and supersonic 
flows that could be manifested as filaments downstream 
of the front. The presence of cells could also increase the 
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area of the front and consequently its global propagation 
rate. If the increase of the flame area is itself accelerative, 
then it could constitute an additional mechanism for the 
observed accelerative expansion of the nebula. 

It is also of interest to investigate whether the detona- 
tion structure suggested here exists in other supernova 
remnants (SNRs), such as those well-observed SNRs in 
different evolution phases. Applications of detonation 
and deflagration theories to describe the chemical or nu- 
clear processes in all kinds of ISM shocks are also merited. 

This work was supported by the startup fund for the 
Center for Combustion Energy at Tsinghua University. 
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